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Correlating Fermi gamma-ray sources with ultra-high 
energy cosmic rays 



o 

(N 



(N 



6 

o3 



> 
oo 
m 

(N 

(N 
O 
O 



N. Mirabal 1 ' 2 * and I. Oya 2 

1 Ramon y Cajal Fellow 

2 Dpto. de Fisica Atomica, Molecular y Nuclear, Universidad Complutense de Madrid, Spain 



ABSTRACT 

The origin of ultra- high energy cosmic rays (UHECRs) is one of the enduring mysteries 
of high-energy astrophysics. To investigate this, we cross-correlate the recently released 
Fermi Large Area Telescope First Source Catalog (1FGL) with the public sample of 
UHECRs made available by the Pierre Auger collaboration. Of the 27 UHECRs in 
the sample, we find 12 events that arrived within 3.1° of Fermi sources. However, we 
find similar or larger number of matches in 63 out of 100 artificial UHECR samples 
constructed using positions randomly drawn from the BATSE 4B catalog of gamma- 
ray bursts (GRBs) collected from 1991 until 1996. Based on our analysis, we find no 
evidence that UHECRs are associated with Fermi sources. We conclude with some 
remarks about the astrophysical origin of cosmic rays. 
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1 INTRODUCTION 

Ultra-high energy cosmic rays (UHECRs) are energetic par- 
ticles (> 10 19 eV) that must originate in the most powerful 
particle accelerators in the Universe. These extreme events 
have fascinated scie ntists from the time of their discovery in 
1962 <|Linslevlll963l). Since then, speculation about their ori- 
gin has nourished JPierpaoli fc Ferrar| [2005; Ghi sellini et al.l 
120081 : ICuesta fc Pradall2009l ). Theoretically, active galactic 
nuclei (AGN) have long been favored as the best candi- 
dates for particle acceleration to these e xtreme energies 
jGinzburg fc Svrovatskiilll964l : lHillaslll984l ). Unfortunately, 
no firm astrophysical association has been established. 

Possibly the most intriguing suggestion of an asso- 
ciation between UHECRs and AGN was put forwa rd by 
the Pierre Auger collaboration II Abraham et al.l [2004) that 
found a possible correlation between the arrival direction 
of 27 events collected in their Southern Observatory and 
the position of nearby ($S 75 Mpc) A GN from the Veron- 
Cetty catalog ( Abrah am et all 120071 ). The correlation re- 
mains but has weakened slightly with the inclusion of addi- 
tional UHECR events (for a grand to tal of 58) collected be - 
tween 2004 January and 2009 March (| Abraham et al.| [2009). 
Nonetheless, questi ons remain about the likelihood of the re- 
ported correlation (|Abbasi et al.|[2008al ). 

It is expected that if AGN are responsible for UHE- 
CRs, these must be preferentially placed within a dis- 
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tance of 100 Mpc (I Abraham et al.ll2007l ). Interactions with 
Cosmic Microwave Background (CMB) photons should 
starve UHECRs arriving from larger dist ances through 
the Greisen-Zatsepin-Kuz min (GZK) effect (jGreisen 1119661 : 
IZatsepin fc Kuzmmlll966T ). In fact, a combination of recent 
observations appear to corroborate the existence of a sup- 
pression of U HECRs above 4 x 10 19 eV consistent with a 
GZK horizon jAbbasi et all l2008bl : [Abraham et al.ll2008bl ). 
It is important to note that recent results from the Pierre 
Auger collaboration point to a transition to heavier compo- 
sition with increasing energies (|Ungerll2007f ). If indeed iron 
nuclei are the dominant component of UHECRs, the large 
deviation angles expected for heavy nuclei would render the 
detection of astrop hysical counterparts ev en the nearest ones 
nearly impossible ([Abraham et alj|2008al ). 

Here, we present a cross-correlation study of UHECRs 
and the recently released Fermi Large Area Telescope First 
Source Cat alog (1FGL) in the 100 MeV to 100 GeV en- 
ergy range (|Abdoet alll2010al ). In contrast with previous 
studies, our cross-correlation analysis is unique in that the 
Fermi catalog comprises a wide range of bona fide particle 
accelerators with gamma-ray production above 100 MeV di- 
rected in the Earth's direction. The Large Area Telescope 
on board Fermi offers a maj or improvemen t in sensitivity 
over previous GeV detectors (|Atwoodl 2007). In its survey 
observation mode, the LAT observes the entire sky every 3 
hours. For individual sources, Fermi provides nearly uniform 
sky coverage down to a photon flux of 4 x 10 -10 cm -2 s _1 
between 1 and 100 GeV, except for sources at low Galactic 
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Figure 1. Distribution of total matches with the 1FGL catalog 
for 100 random sample drawn from the BATSE 4B catalog. The 
vertical line indicates the number of matches between UHECRs 
detected by the Pierre Auger Observatory and the 1FGL catalog. 



latitude (\b\ ^ 10° ) where the diffuse emission dominates 
l|Abdo et al.ll2010ah . Throughout the paper, we assume an 
H =71 km s" 1 Mpc~\ fl m = 0.27, Qa = 0.73 cosmological 
model. 



2 CROSS-CORRELATION OF Pierre Auger 
UHECRS WITH Fermi SOURCES 

As of 2007 November, the public catalog of UHECRs 
released by the Pierre Auger collaboration contains 27 
events with energies ab ove 55 EeV collected at their site 
in Malargiie, Argentina (| Abraham et al.ll2004h . The Fermi 
LAT 1FGL catalog consists of 1451 sou rces characterized i n 
the 100 MeV-100 GeV energy range (|Abdo et all |2010ah . 
The data were obtained in an all-sky scanning mode during 
2008 August-2009 July and represents the most extensive 
map of the gamma-ray sky 100 MeV) ever obtained. The 
entire catalog includes 689 blazars, two starburst galaxies, 
two radio galaxies, 56 pu lsars, 50 supernova remnants, and 
630 unidentified sources (lAbdo et al.ll2010al ). 

Since we are interested in testing for possible correla- 
tions of UHECRs with various classes of gamma-ray emit- 
ters in the MeV-GeV energy range without any a priori as- 
sumption, we take advantage of the complete 1FGL catalog 
without any redshift or type discrimination. 



To test for a possible cross-correlation between UHE- 
CRs and Fermi sources in the 1FGL catalog, we checked 
whether the individual UHECR positions in the Pierre 
Auger sample are clearly contained within the error circle of 
individual sources listed in the 1FGL catalog. Specifically, 
we count correlations whenever an UHECR event is within a 
circle of 3.1° radius around a particular Fermi source. This 
radius is in line with the value predicted by conventional 
models of cosmic ray trajectories t hat consider the full e f- 
fect of the Galactic magnetic field (| Abraham et al.ll2008al ). 
In order to avoid "multiple" counts, we only consider one 
match per individual UHECR event. 

Within the 27 UHECR events, we find 12 matches with 
1FGL sources. To examine the likelihood of such correla- 
tion, we used the BATSE 4B gamma-ray burst (GRB) cat- 
alog (|Paciesas et alj|l999h in place of a random generator 
of isotropic sky positions that allows us to generate artifi- 
cial samples of simulated UHECRs. The BATSE 4B catalog 
consists of 1637 GRB positions localized by the BATSE in- 
strument on board the Compton Gamma Ray Observatory 
(CGRO) in the period between 1991 April 19 and 1996 Au- 
gust 29. For our work, we restricted our analysis to events 
that would be accessible from the Pierre Auger southern site 
at a declination < 24.8° (| Abraham et al-lliooil ). In addition, 
we assume that the fraction of exp osure is the sam e accros s 
the declination range as stated in lAbraham et al.1 (|2008al ). 
After the proper cuts were applied, we drew 100 random 
sets of 27 events from the resulting "southern" BATSE 4B 
catalog to match the UHECR sample. We next proceeded to 
correlate each of the 100 random sets with the 1451 Fermi 
sources in the same manner as with the original Pierre Auger 
dataset. 

Figure Q] shows the distribution of matches between the 
artificial samples of UHECRs constructed from the BATSE 
4B catalog and the 1FGL catalog. In particular, 63% of the 
artificial samples have 12 or more matches consistent with 
Fermi positions. Allowing for larger circles (6° and 8° radii) 
only increases the number of false positives when compared 
with the BATSE sample to 71% and 73% respectively. If 
the position of Pierre Auger events were strongly correlated 
with Fermi sources, we would expect a greater number of co- 
incidences when compared with the outcomes of randomly- 
generated artificial sets of UHECR events drawn from the 
BATSE 4B GRB catalog. We find no evidence for such an 
association and subsequently cannot claim a positive cross- 
correlation between UHECRs and the 1FGL catalog. 



3 DISCUSSION AND CONCLUSIONS 

In summary, we find no cross-correlation between UHE- 
CRs and 1FGL sources that cannot be reproduced by 
ch ance alignment. This d iffers from the findings reported 
by lAbraham et~afl J2009I). using a different AGN sample 
(see also Abraham et alj 20071 ). Examining the matches be- 
tween the Pierre Auger events and the 1FGL catalog in 
closer detail, we notice that the sample includes one pulsar, 
two unidentified sources, and eight AGN. Four of the AGN 
have measured well beyond the GZK horizon with redshifts 
as high as z = 1.843. The lowest redshifts correspond to 
Centaurus A (z = 0.002) and NGC 4945 (z = 0.002), as 
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we shall discuss later. The remaining two AGN associations 
have no spectroscopic redshifts. 

One could argue that the aforementioned AGN without 
spectroscopic redshifts could potentially form a parent pop- 
ulation for UHECRs. If it was, these AGN would most likel y 
lie at distances closer than 100 Mpc (I Abraham et alj|2008al ). 
Optical imaging studies of BL Lac reveal that the host galax- 
ies of nearby AGN with jets pointed in our direction tend to 
cluster around a mean absolu te magnitude in the R band 
equiv alent to Mr = —22.8 (|Sbarufatti. Treves fc Falomol 
120051 ). For distances closer than 100 Mpc, this would cor- 
respond to an apparent i?-band magnitude R = 12.2 or 
brighter. A quick analysis of optical Digitized Sky Survey 
images centered on the positions of the two matched AGN 
without spectroscopic redshifts reveals no optical sources 
brighter than R = 12.2 at the derived AGN positions 
l|Abdo et al.ll2010bl ). Therefore, it is unlikely that any of the 
paired UHECRs/ AGN are located within the hypothesized 
GZK horizon at 100 Mpc. 

Interestingly, we do recover the two UHECRs pre- 
viou sly associated with the radio galaxy Centaurus 
A ijAbraham et alJ [2008a]) and one UHECR consis- 
tent with the position of th e starburst/Seyfert 2 NGC 
4945 (jMoskalenko et all [20091 ). Both Centaurus A (1FGL 
J1325. 6-4300) and NGC 4945 (1FGL J1305.4-4928) are de- 
tected by Fermi. However, NGC 4945 only stands out as one 
of the two Seyfert 2 galaxies in the 1FGL catalog rather than 
for its qualifications as a UHECR accelerator ()Abdo et al.l 
2010b). On the other hand, radio galaxies such as Centau- 
rus A could potentially account f or the acceleration o f some 
UHECRs to energies > 55 EeV (jDermer et al.ll2009l ). As a 
result, it is important to continue to explore the connection 
(if any) between UHECRs and Centaurus A. However, it is 
troubling that 24 UHECR events detected by Pierre Auger 
remain without an apparent Fermi gamma-ray counterpart 
within 100 Mpc (we dub these "orphan" UHECRs). 

In fact, there seems to be an apparent dearth of nearby 
AGN detected in gamma rays at distances less than 100 
Mpc (z ^ 0.025). To see this more clearly, Figure [2] shows 
the distr ibution for Fermi A GN in the redshift range z = 
- 0.1 ijAbdo et al.l l2010bl ). In total, there are only five 
objects with z ^ 0.025 that would be accessible to the Pierre 
Auger southern site namely NGC 253 (z = 0.001), NGC 
4945 (z = 0.002), Centaurus A (z = 0.002), M 87 (z = 
0.004), and ESO 323-G77 (z = 0.015). NGC 4945 (being 
generous) and Centaurus A could potentially account for 3 
of the 27 UHECR events as discussed above. However, we 
remark that there are no additional identified Fermi sources 
in the 100 MeV-100 GeV energy band to account for the bulk 
of UHECRs seen by Pierre Auger so far. 

It is possible that other AGN remain unidentified in 
the 1FGL catalog, as implied by the observed north-south 
anisotropy of the associated sources, which might reflect the 
incompleteness of t he existing AGN catalogs in other bands 
l|Abdo et alJl2010bl ). However, the remainder is expected to 
lie in the lower end of the gamma-ray flux distribution, and 
therefore in principle not represent the most efficient particle 
accelerators. In addition, some AGN are known to be highly 
variable, therefore we cannot fully discard that some sources 
have escaped Fermi detection and possible correlation with 
UHECRs because of a low state. 

Accordingly, one must admit one or all of the following 




AGN redshift 

Figure 2. Redshift distribution of Fermi AGN between z = 
and z = 0.1. Northern (deel > 24. g 3 ) and southern (deel < 24. f?) 
samples are shown in blue (dashed line) and grey (solid line) 
respectively. Only five identified Fermi sources have z < 0.025 
and are accessible to the Pierre Auger southern site. The vertical 
line marks a distance of 100 Mpc. 

possibilities: 1) if UHECRs are truly accelerated within the 
current sample of detected Fermi sources, their trajectories 
must experience deviations greater than 3.1° from their ac- 
tual astrophysical origin, making the identification of UHE- 
CRs nearly intractable. In particular, the absence of astro- 
physical counterparts could be well justified if UHECRs are 
dominated by heavy nuclei including iron since the devi- 
ation scales linearly the atomic number Z, 2) there could 
potentially exist an undetected p opulation of near by AGN 
- the so called proton blazars? (|Mannheiml Il993l ) - that 
only emit in a very narrow band above 100 GeV to be dis- 
covered by future VHE experiments such as the Cherenkov 
Telescope Array (CTA) or the Advanced Gamma Imaging 
System (AGIS), and 3) the absence of obvious particle accel- 
erators in the MeV-GeV energy band leaves some room for 
the exploration of additional forms of particle acceleration 
and even exotic possibilities. 

Trying to narrow down future observational directions, 
it is important to search for possible timing/directional cor- 
relations of the arrival direction of cosmic ray with AGN 
flares. Such exercise currently lacks the dedicated sky mon- 
itoring that could potentially identify the majority of AGN 
flares in real time. It would also require the rapid release of 
UHECR detections (including position) in real time. The ex- 
isting fleet of MeV-TeV telescopes including AGILE, Fermi 
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H.E.S.S., MAGIC, and VERITAS could help trace strong 
flares in the gamma-ray band. However, groundbreaking 
progress into the gamma-ray variability domain might have 
to wait for the next generation of VHE experiments. In the 
theoretical front, it is critical to improve existing models 
of cosmic ray propagation. In particular, based on current 
cross-correlation studies, the initial expectation that cosmic 
ray trajectories at energies > 10 19 eV should be fairly rigid 
does not appear to be so obvious. 
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